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ABSTRACT 

We attempt to c onfirm bright non-L TE emission from the exoplanet HD 189733b at 3.25 /im, as 



recently reported by [Swain et al. (2010 ) based on observations at low spectral resolving power (X/SX 
30). Non-LTE emission lines from gas in an exoplanet atmosphere will not be significantly broadened 
by collisions, so the measured emission intensity per resolution element must be substantially brighter 
when observed at high spectral resolving power. We observed the planet before, during, and after 
a secondary eclipse event at a resolving power X/SX = 27,000 using the NIRSPEC spectrometer on 
the Keck II telescope. Our spectra cover a spectral window near the peak found by Swain et al., and 
we compare emission cases that could account for the magnitude and wavelength dependence of the 
Swain et al. result with our final spectral residuals. To model the expected line emission, we use a 
general non-equilibrium formulation to synthesize emission features from all plausible molecules that 
emit in this spectral region. In every case, we detect no line emission to a high degree of confidence. 
After considering possible explanations for the Swain et al. results and the disparity with our own 
data, we conclude that an astrophysical source for the putative non-LTE emission is unlikely. We note 
that the wavelength dependence of the signal seen by Swain et al. closely matches the 2^2 band of 
water vapor at 300K, and we suggest that an imperfect correction for telluric water is the source of 
the feature claimed by Swain et al. 

Subject headings: astrobiology - infrared: planetary systems - planets and satellites: individual 
(HD189733b) - radiative transfer - techniques: spectroscopic 



1. INTRODUCTION 

Characterization of exoplanet atmospheres is a highly 
active field of research, especially using transit and 
eclipse techniques. Space-borne detections of transiting 
exoplanet atmospher es have become commonplace, both 



Deming 



in the infrared (e.g., |Charbonneau et al.||2005 
et al.||2005| |Harrington et al.||2007| |C harbon neau et al. 
2008| [Swain et al.||2008; Knutson et al.||20091 |Pont et a 



20081 |2009[ IGrillmair et al.||2008| |Fressin et al.||2010 
Gillon et ai l 1 2010 i and optical spectral regions (e.g., 
Charbonneau et al.|2002"||Snellen et al.|2009[|2010b| ). The 

space- borne data have been extensively exploited to in- 
vestigate the temperature str ucture and dynamics of hot 



exoplanet atmospheres (e.g., Burrows et aL]p007, 2008; 



Fortney et aT1[2008l I Showman et 
2010 sec Seager & Deming 2010 



al. ||2009||Fortney et al.| 



for a recent review) 



Near-intrared photometry of the hottest exoplanets 
is now possible using ground-based observatories (e.g., 



Gillon et al.|2009] |Alonso et alT20To] |de Mooij fc Snellen| 



2009p , and spectroscopy of strong atomi c lines has also 



been successful from the g round (e.g., Redfield et al. 



2008 Snellen et al. 2008). Since a new generation 



ot extremely large ground-based telescopes is now be- 
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ing planned (Hook 2009), ground-based spectroscopy of 



molecular features m exoplanet atmospheres may even- 
tually become possible, and would provide exceptional 
diagnostic power. So far, ground-based attempts to 
detect near-IR molecular features in exoplanet atmo- 
spheres with a range of spectral resolving powers have 



for the most part been unsuccessful (e.g., Wiedemann 



et al.|2001HBrown et al. 2002; Ri chardson et al.|2^03aTp] 



Barnes et al.||2007[ |2008j]2010|) , though high-resolution 
spectroscopy has recently produce d a successful detec - 
tion of CO at K-band wavelengths ( Snellen et al.|2010a ) 



Additionally, a startling ground- based detection of 
molecular spectral features from the atmosphere of the 
exoplanet HD 189733b at moderate spectral resolution 
was announced by Swain et al. (2010; hereafter S10). 
S10 analyzed observations of HD 189733 taken with the 
SpeX instrument on the Infrared Telescope Facility, fo- 
cusing on two wavelength regions: 2.0 — 2.4 /im (K band) 
and 3.1 — 4.1 /im (L band). Their results show a rising 
spectral slope at 2.2 /im, which they attribute to CO2 
absorption, and a bright peak at 3.25 /im. The shape 
of the spectrum at K-band wavelengths matched previ- 
ous results at the se wavelengths tak en with the Hubble 



Space Telescope (Swain et al. 2009), which the authors 
construed as support tor their ground-based detection in 
the L band. 

S10 interpret the flux peak in the L-band portion of 
their data as evidence of bright non-LTE (NLTE) emis- 
sion from radiatively-excited methane. The wavelength 
region for the claimed emission roughly corresponds with 
the R-branch region of CH 4 (3.1 — 3.34 /im), but no cor- 
responding feature is seen in the region of the methane 
P-branch (3.35—3.5 /im); S10 note this disparity but offer 
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no theoretical explanation for it. The L-band emission 
claimed by S10 is also remarkable for having high inten- 
sity at low spectral resolution. The wavelength extent of 
each spectral bin analyzed by S10 is 0.1 /im, (A/<5A ~ 30), 
and they find a planet-to-star contrast ratio of approxi- 
mately 0.9% in their brightest bin at 3.25 /im. This re- 
sult is even more surprising considering that NLTE emis- 
sion lines arc unlikely to be collisionally broadened, since 
the pressures needed to achieve significant line broaden- 
ing would produce high collision rates that would drive 
the level populations to thermal equilibrium. Therefore 
any NLTE emission lines will have maximum intrinsic 
widths comparable to the Doppler velocities of the emit- 
ting species, and when observed at low resolving power 
the apparent flux densities (ergs s _1 cm -2 Hz -1 ) of emis- 
sion lines will be greatly reduced because their intrin- 
sic line widths are unresolved. This motivated our at- 
tempt to confirm the S10 results using data with a much 
higher spectral resolving pow er, obtained with th e Near- 
IR Spectrograph (NIRSPEC; |McLean et al.|1998| on the 
10-m Keck II telescope . The peak intensities of NLTE 
emission lines will be much greater when they are ob- 
served at high spectral resolving power, and SlO's claim 
of bright emission at low spectral resolution implies that 
observations at high spectral resolution and sensitivity 
should easily produce a detection. 

A description of our search for absorption signatures of 
the exoplanet in our data set will be presented in a future 
manuscript; the expected line depths based on standard 
atmospheric modeling are approximately .01 — 0.1% of 
the stellar continuum ( Burrows et al.|[2008| , and a com- 
prehensive analysis at that level is beyond the scope of 
this study. In this paper we limit our analysis and discus- 
sion to whether or not we can confirm the bright emis- 
sion signal announced by S10. We present our observ- 
ing and data analysis procedure in Section £j2j predict 
the expected emission features using excitation models of 
molecular emission at different rotational temperatures 
and compare the expected signal to our results in Sec- 
tion Ej3j and conclude with a discussion of the potential 
explanations for the differences between our results and 
those of S10 in Section |4] 

2. OBSERVATIONS AND DATA ANALYSIS 

We acquired spectra at A/ A A « 27,000 with NIRSPEC 
without AO on UT July 13, 2009 using the KL filter with 
a setting covering portions of the wavelength range be- 
tween 3.27 and 4.0 /im. Observing conditions were opti- 
mal, with low water vapor and clear skies. We observed 
HD 189733 before, during and after a secondary eclipse 
of the planet, for a total integration time of 100 min- 
utes between UT 10:00 and UT 14:00. A bright B-type 
comparison star (HR8634) was also observed immedi- 
ately after the science target. We nodded the telescope 
12 arcsec in an ABBA sequence, with 60-second integra- 
tions per beam for both stars. In total, we obtained 48 
cchcllc spectra of HD 189733 during eclipse, 52 spectra 
out of eclipse, and 40 spectra of HR8634. 

Since our spectra were acquired prior to the pub- 
lication of the S10 results, the wavelength ranges of 
our echelle grating orders are not identical to the S10 
work. Fortunately, our data include a spectral range from 
3.27 — 3.31 fim that overlaps the brightest bin in the emis- 
sion feature claimed by S10. The wavelength structure of 



molecular bands is known unequivocally from quantum 
mechanics, with only the level populations affecting the 
intensity of lines in different spectral channels. Our mod- 
eling (see Section ^3| indicates that molecules emitting 
strongly in the S10~3.25/im bin must also emit signifi- 
cantly in our 3.27 — 3.31 fim region, and we test the S10 
results on that basis. 

We utilized custom data reduction algorithms, previ- 
ously used to detect new molecular e mission features 
from warm gas in circumstellar disks (Mandell et al. 



20081, to extract and process spectra lor each echelle 
order in each ABBA set. We reduced the initial 2D 
spectral-spatial images to ID spectra after first correct- 
ing for the slope of the beam due to cross-dispersion and 
subtracting A- and B-beam images to remove the contri- 
bution from telluric radiance. We identified bad pixels 
and cosmic ray hits in each raw pixel column by compar- 
ing the beam profile to an average beam profile for nearby 
columns, allowing us to identify and replace single-pixel 
events without removing any enhancements due to emis- 
sion or absorption features. 

We corrected for changing airmass and telluric at- 
mospheric conditions to high accuracy using a two- 
step process: 1) fitting the data for both the science 
star (HD 189733) and the comparison star (HR8634) 
with terrestrial spectral transmittance m odels synthe- 
sized with the LBLRTM atmospheric code ( |Clough et al.| 
2005 1 and subtracting the models to obtain spectral 
residuals for each star, and 2) differencing the residuals 
of the two stars. The subtraction of the telluric model 
compensates for the effects of changing airmass and at- 
mospheric variability, and the differencing of the residu- 
als removes remnant fringes and other instrumental ar- 
tifacts, as well as minor errors in the telluric model such 
as imprecise pressure broadening, weak features missing 
from the line list, and inaccurate isotopic ratios. For 
our LBLRTM models we utilized a standard tropical 
temperature profile and line parameters from the HI- 
TRAN2008 molecular database with updates from 2009 
( Rothman et al. 2009 ), and we fitted for the abundance of 
three key atmospheric constituents (H2O, CH4, and O3) 
and a scaling factor for the temperature of the tropo- 
sphere. Atmospheric models were fitted for wavelength 
sub-sections of each AB set, with the fully- resolved model 
convolved with the local instrumental resolving power 
derived for each sub-section to compensate for variabil- 
ity in the effective spectral resolving power due to the 
position of the spectrum on the detector. Additionally, 
using a telluric model provided us with an extremely 
well-calibrated wavelength solution for each sub-section 
(AA/A ~ 10 -6 ). This process achieved results corre- 
sponding to S /N ~ 300 on the original stellar continuum 
for each AB set, and an rms noise only slightly larger (by 
20%) than that expected from the photon statistics (see 
Figure [l]) . 

Stellar absorption features were removed by using the 
results for the in-eclipse data as a stellar template, with 
each AB set shifted to correct for changes in the heliocen- 
tric velocity and the stellar radial velocity variation due 
to the influence of the planet. We then subtracted this 
stellar template from each AB set in the out-of-eclipse 
data after shifting for changes in the radial velocity for 
that time interval, leaving final residuals for each out-of- 
eclipse AB set with only the signal from the planet (see 
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Fig. 1. — Demonstration of our data analysis procedure. Synthetic terrestrial atmospheric models are fitted separately to both the 
science and comparison stars and subtracted (1); then the two sets of residuals are differenced to remove second-order instrumental or 
atmospheric features (2). Remaining fringing is removed using a Fourier filter (3). The S/N in the final residuals is ~ 1000. Stellar features 
are removed by observing the star during the eclipse of the planet and using the data to create a template of the stellar spectrum which is 
then subtracted from the out-of-eclipse data. 

Figure [2]). 

We combined all the out-of-eclipse residuals after shift- 
ing each set to correct for the radial velocity of the planet 
around the star, wit h velocity shifts calcula ted based on 
the ephemeris from IKnutson et al. (2009). The radial 
velocity of the HD 189733 system relative to the tele- 
scope (accounting for Earth's orbit and rotation) was 
-8 km/s, but over our observing window the exoplanet's 
orbit led to velocity shifts over 60 km/s. Because the 
data must be shifted in velocity space to account for 
the shift of the planet's spectrum, sections of the spec- 
trum with high transmittance may be shifted to line 
up with spectral regions of low transmittance; this ben- 
efits us by moderating extinction of lines with strong 
telluric counterparts. However, to avoid degrading our 
detectable signal by combining low-transmittance data 
with high-transmittance data, we weighted the value for 
each shifted spectral channel by the stochastic-noise SNR 
(based on detected photons and transmittance at the pre- 
shifted wavelength). 

Finally, we applied a high-pass filter to the residuals 
for each set to remove variations in the continuum due 
to movement of the star on the detector over the course of 
the night by creating a smoothed version of the spectrum 
using a boxcar average with a window sufficiently large 
that any narrow features would be significantly reduced 
(~ 20 resolving elements), and then subtracted the fil- 
tered spectrum to achieve a final rms scatter of ^0.001. 
We evaluated the loss of narrow features due to the filter- 
ing by injecting synthetic emission lines and determined 
that the high pass filter preserved 96% of their ampli- 
tudes. Flux calibration was performed by normalizing 
the observed continuum flux to t he predicted flux b ased 

[ 2003]) and 
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Fig. 2. — Sample spectral region for a single AB data set for 
HD 189733. A telluric model (green) is fitted to the observed spec- 
trum and subtracted to produce an initial set of residuals. We then 
subtract residuals from a similarly processed calibration star (blue) 
to remove instrumental errors and systematic deviations from the 
telluric model. We then bin the in-eclipse data to create a stellar 
template (red) which we subtract from the out-of-eclipse data to 
achieve our final residuals. Low-amplitude deviations remain in 
the final residuals, primarily due to changes in the instrumental 
resolving power, but these features are eliminated when multiple 
sets are combined due to the changing radial velocity shift. 

3. A SEARCH FOR MOLECULAR EMISSION FEATURES 
AT HIGH RESOLUTION 

In order to understand the origin of the excess emis- 
sion signal detected by S10 between 3.0 and 3.5 /im, we 
must first understand the potential contributing emit- 
ting species and their spectral characteristics. S10 sug- 
gest that methane emission at 3.3 /jm provides the best 
fit to their data, but there are clearly other species with 
transitions in this spectral region that may also be suf- 
ficiently abundant to produce detectable emission in the 
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planet's atmosphere. In order to constrain the potential 
molecular species that could be producing emission, we 
generated models of line emission for a range of possible 
molecular species and compared the resulting predictions 
to our data. The goal of these models is not to explain 
the S10 results through a prediction of emission from a 
specific mechanism; rather, we only want to mimic the 
location and magnitude of any emission lines that might 
arise from a variety of molecules assuming the most basic 
physical requirements, under a range of excitation con- 
ditions, and then scale this emission based on the S10 
results. 

The dominant gas-phase molecular features in the 
L band are ro-vibrational spectral lines, in which a 
molecule transitions from a specific higher-energy rota- 
tional and vibrational configuration to a lower-energy 
rotational and vibrational state. The population of 
molecules in any specific rotational and vibrational state 
is determined by the transition rates between vibrational 
energy states as well as the transition rates between rota- 
tional energy states. The fundamental ratios between the 
strengths of different transitions for a gas in local ther- 
modynamic equilibrium (LTE) are set by the Boltzmann 
distribution and are based on the temperature of the gas 
and the probability for each individual transition. Any 
NLTE mechanisms such as fluorescence (radiative pump- 
ing by incident photons) or charged-particle interactions 
will increase the excitation rates and subsequent emission 
intensities; the level populations are therefore a function 
of both the underlying equilibrium distribution as well as 
any NLTE excitation processes. In particular, fluorescent 
excitation of irradiated gas may result in an increased 
strength for higher-energy vibrational bands whose up- 
per state populations are pumped by high-energy pho- 
tons from a nearby radiation source (e.g., the central 
star). 

For comparisons with the data, we do not want to as- 
sume a specific physical excitation mechanism a priori 
(we discuss the potential for different excitation mech- 
anisms in more detail in Section Sj4j), but we require 
that the basic rotational and vibrational structure of 
each molecule's spectrum is maintained. Since transi- 
tions between rotational energy states require much less 
energy than vibrational transitions, we first examined 
whether rotational transitions for energy states in the L 
band would be dominated by LTE or NLTE excitation. 
Using excitation rates from the HITR AN2008 molecu- 



lar database with updates from 2009 (Rothman et al 



2009 1 , and collisional rate coefficients for rotational tran- 



sitions of our candida te molecules with molecular hy- 



drogen (10 11 cm 3 /s, Faure et al. 2008 Faure fc Jos 



selin|2 008 ) and free electrons (10 b — 10 9 cm J /s, Faure 



et al. 2004 1 , rotational thermalization of the ground vi 



brational states is expected to be achieved at pressures 
greater than 0.01 millibars. Only at very high altitudes, 
where total column densities are negligible and many 



molec ular species would be photodissociated ( Liang et al. 
|2003[ ), does rotational equilibrium fail. Therefore, the 
rotational structure for any ro-vibrational emission from 
the exoplanet will broadly resemble a Boltzmann distri- 
bution for the rotational population, thermalized to the 
local kinetic temperature, and any additional NLTE exci- 
tation will only affect the intensity of different vibrational 



energy states. 

We therefore calculate relative line strengths for ro- 
vibrational transitions for a candidate molecule using two 
bracketing cases: 1) we adopt the simplest case of a com- 
mon rotational and vibrational excitation temperature 
(i.e., all vibrational bands are collisionally excited) from 
an optically thin gas, and 2) we postulate that emis- 
sion from a single radiatively excited vibrational band 
dominates the spectrum. This formulation is quite gen- 
eral, and allows us to test the results of exciting a spe- 
cific vibrational band without prior assumptions about 
the relative vibrational level populations. We first cal- 
culated rotational and vibrational levels populated in 
proportion to a single characteristic excitation temper- 
ature that we varied over a wide range (from 100K to 
10,000K, stepping by 0.2 in logio space); examples are 
compared with the S10 results in Figure [3j We then cal- 
culated models with only a single vibrational level popu- 
lated. For each of these cases, we calculate the emission 
spectrum o f each g i ven m olecule, and divide that spec- 
trum by a Kurucz (1993) model stellar spectrum that 
matches the star (T e// = 5000K, log(g)= 4.5, Z= 0.0). 
We adjust the absolute level populations of the emit- 
ting molecule so that the resultant contrast (planet/star) 
equals the observed S10 value at 3.25 /im, and we use a 
X 2 fit to the five data points between 3.05 and 3.45 fim 
from S10 to determine the best-fit rotational temper- 
ature for both the single-temperature models and the 
single- vibrational-level models. We then examine the im- 
plied molecular spectrum in our adjacent wavelengths to 
determine whether our spectra are consistent with that 
molecule at that excitation temperature. 

We determined that the species with transitions in this 
region, and the potential for a reasonable number of emit- 
ting molecules, are H2O, CH4, NH3, and OH; modeling 
quickly shows that NH 3 peaks below 3.25 /im for all tem- 
peratures and vibrational levels. In Figure [3j we plot 
the best-fit single-temperature models for the other three 
molecules (475K for CH 4 , 300K for H 2 and 500K for 
OH) , as well as two additional single-temperature models 
covering a large range in excitation temperature. We also 
plot the best-fit single-vibrational-band model for each 
molecule. In the lower panels, the models are binned in 
the same spectral channels as S10 and scaled to match 
their data point at 3.25 /im (their bin overlapping our 
measured wavelength region). 

It is clear from our modeling that none of the single- 
temperature models for CH4 and OH fit the S10 results. 
The single-vibrational-level models allow a reasonably 
good fit, though each model has at least one major devi- 
ation from the S10 data. Additionally, the scaling factors 
required to fit the CH4 and OH single-vibrational-level 
models to the data (10 2 and 10 9 respectively) require a 
level of radiative pumping that would be unattainable 
based on the stellar flux, as we explain belo w when con- 
sidering possible fluorescence (Section [4.1). There may 
be an intermediate model that fits the S10 data with even 
more precision, with a number of vibrational bands at a 
range of intensities, but we are mostly interested in the 
prediction for emission features in our spectrum and the 
bracketing cases considered here provide the full range 
of possible features; a full NLTE model is beyond the 
scope of this paper. Interestingly, the models for H2O 
are able to fit the data very well (\ 2 ~ 1) using low- 
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Fig. 3. — Excitation models for the three candidate molecular species that wc compare to our data to search for expected emission 
features. For each molecule, we calculate three models using a single excitation temperature for the Boltzmann distribution of both the 
rotational levels and the vibrational levels. In the fourth panel on the right, we plot the best-fit model using only a single vibrational 
level, in order to test th e potential effects o f strong radiative pumping of an upper vibrational level. We scale each model to the 3.25 fim 
planet-star contrast fromlSwain et al.|(|2010|, and overplot their results in red. H2O clearly provides the best fit . The single-yibratio nal-level 
models for CH4 and OH also ht reasonably well, but the scaling required to reproduce the flux reported by Swain et al.| pjip) requires 
unphysical pumping efficiencies. 



temperature, single-excitation-temperature models (the 
single- vibrational-level model fits as well, due to the dom- 
inance of the single 2^2 vibrational band in the spectral 
region) . We disc uss the potential implications of this 
result in Section §4.4| 

In order to run a full grid of models for a range of tem- 
peratures and vibrational states for each molecule over 
a wide wavelength range, it was necessary to utilize a 
list of transitions that was both relatively complete and 
manageable in size. We primarily used data from the 
HITRAN2008 database for the full grid of models; how- 



ever, we then compared our results to models computed 
with more complete line lis ts - data from the CT LiQTitan 
line database for methane ( Boudoir et al.|20"To]) and data 
from the HITEMP20I0 database ( jRothman et al.||2010 | 
for water and OH. Differences between the low-resolution 
models based on HITRAN and the more complete models 
were virtually non-existent at lower temperatures, and 
while the flux for some increased up to 20% for the high- 
est temperatures, these differences are insignificant due 
to the already-poor fit between higher-temperature mod- 
els and the S10 results as illustrated in Figure [3] and the 
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HD 1 89733, NIRSPEC/Keck, Order 23, July 2009 
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Fig. 4. — Results from our data reduction, with candidate molecular models using the CH4@Titan and HITEMP2010 line lists overplottcd. 
The upper black trace is the original data, and the bottom black trace are our residuals after removing telluric and stellar features, shifting 
the data to the correct radial velocity and combining it, and then removing large-scale gradients in the continuum with a high-pass filter 
(as described in Jf2J) . Models for CH4 (top) and H2O (bottom) with different excitation temperatures are plotted above. We detect no 
emission at any ofthe expected positions, with upper limits more than 10 times below the total intensity of the expected emission for each 
model. 

fact that our upper limits are significantly smaller than 
the expected signal. 

In Figure |4j we plot the scaled models based on the 
more complete data sets over our spectrum, which cov- 
ers both high- and low-temperature methane and water 
transitions. The models are multiplied by the telluric 
transmittance model (i.e., attenuated by extinction from 
the Earth's atmosphere) and filtered using the same high- 
pass filter used to remove the continuum in the data. 
None of the candidate molecular emission features that 
could account for the S10 result can be detected in our 
data. Indeed, in every instance they can be ruled out 
unequivocally; in most cases the emission would be so 
bright that the features would be visible in the raw data. 
We also performed the same data reduction procedure 



on three additional NIRSPEC orders, spanning wave- 
length ranges of 3.41 — 3.46 /im, 3.57 — 3.63 /im, and 
3.75 — 3.81 /im; no emission was detected to the same 
sensitivity limits in any ofthe orders. The standard devi- 
ation of our observed residuals (lower trace in Figure[4| is 
0.0011 at our observed resolution, and is even lower when 
averaged over a line width. To see how firmly the vari- 
ous cases are rejected, consider the spectrum of methane 
excited at 1200K. That spectrum requires a 3.5% emis- 
sion feature near 3.27 fim that would be detected at more 
than 30er significance in our data, but is not seen. The 
other molecular cases illustrated in Figure [I] are rejected 
at even higher levels of significance. Moreover, as stated 
above, we have not restricted our search to merely the 
cases that are illustrated in Figure [4] In no instance can 
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we find a spectrum that accounts for the S10 results but 
is also allowed by our data. We can only reconcile our 
data with the S10 result if we allow only the rotational 
transitions not covered by our spectral windows to con- 
tain flux, which we consider impossible given the range of 
transitions with similar excitation energies between 3.0 
and 3.5 fxm as well as the high collisional cross-sections 
for changes in the rotational quantum number of water 
and m ethane ( Faure et al.| |2004 2008 Faure & Josselin 
20081. 



4. DISCUSSION 

In this section we discuss possible explanations for the 
emission feature detected by S10, and evaluate whether 
these explanations are consistent with known physical 
emission mechanisms as well as the observational results 
from this study and previous studies of the HD189733 
system. We then address the possibility that the emission 
feature is not astrophysical in nature, but rather a data 
reduction artifact caused by incomplete removal of the 
signature of telluric water vapor. 

4.1. Fluorescent Pumping by the Stellar Continuum 

The most well-understood NLTE process that can pro- 
duce NIR line emission is fluorescence excited by ra- 
diation from the stellar continuum, either through di- 
rect pumping by incident radiation at resonant frequen- 
cies (resonant fluorescence) or through radiative cascades 
from upper energy levels (non-resonant fluorescence). 
In our own Solar System, methane fluorescence is well 
known, for example from the at mospheres of Jupite r 
prossart et~aT| |1999| ) and Titan flBrown et aT||2006[ ), 
and S10 claim methane fluorescence as the most likely 
source for their signal. 

Both resonant and non-resonant fluorescence pumped 
by photons from the stellar continuum can be rejected as 
the source of the S10 feature on both observational and 
theoretical grounds. First, any emission would be present 
at a similar level in both 2007 and 2009, and we con- 
clusively rule out the presence of bright narrow emission 
features that would be present in our data. Furthermore, 
consideration of the energetics required to produce fluo- 
rescent emission shows that the emission feature claimed 
by S10 is so bright that it is impossible to attribute to 
fluorescence based on pumping by stellar radiation. 

For resonant fluorescence, any radiation scattered back 
to the observer by the planetary atmosphere cannot be 
more intense than the impinging stellar radiation; this 
provides an "optically thick" limit for an absorbing tran- 
sition. The radii of the planet and star, and other system 
parameters, a re well known from high-precision transit 
observations ( Winn et al.||2007[ ). Therefore, we know 
that the solid angle subtended by the planet is 2.1% of 
the stellar solid angle as viewed from Earth. S10 find 
a contrast (planet flux divided by stellar flux) equal to 
0.9% in their brightest bin at 3.25 /im. If we start with 
the assumption that their feature is at least marginally 
resolved (i.e. their measured peak flux is not diluted over 
several resolving elements) , that implies that the specific 
intensity (ergs cm -2 sec -1 sr _1 Hz -1 ) of the planet's 
emission at the top of its atmosphere is 43% of the spe- 
cific intensity of stellar radiation at the top of the stellar 
atmosphere. However, since the planet orbits at 8.4 stel- 
lar radii, the intensity of the stellar radiation field im- 



pinging on the planet is geometrically diluted by a factor 
of 8.4 2 = 70. Hence, in order to produce the brightness 
of the S10 feature, the specific intensity of the emission 
from the planet would have to be 30 times greater than 
the impinging stellar radiation at 3.25 /im to produce the 
peak flux measured at low resolution. 

If we assume that any emission is actually not resolved 
(as would be expected for NLTE emission from a low- 
density medium), the discrepancy becomes even more 
pronounced. For example, consider our molecular exci- 
tation model for methane at a 1200K rotational temper- 
ature, as described above. This case produces a spec- 
trally resolved emission line at 3.27/im that would peak 
at 3.5% of the stellar continuum in our spectrum, but 
without extinction by the terrestrial counterpart would 
actually peak at 13.5%. Adjusting for the ratio of solid 
angles, the specific intensity in the core of this line is 
6.5 times the stellar intensity at the top of the stellar 
atmosphere, at the same wavelength. Since the stellar 
radiation field is diluted by a factor of 70 at the orbital 
distance of the planet, this specific planetary line would 
be 450 times brighter than the impinging stellar radiation 
field at that same wavelength. Resonant fluorescence is 
therefore incapable of producing the apparent flux ob- 
served by S10, whether the measured emission is fully 
resolved by S10 or not. 

Although resonant fluoresence alone is not a feasible 
explanation for the S10 emission feature, we consider 
whether other fluorescence mechanisms could contribute 
significantly to the emission reported by S10. For exam- 
ple, stellar photons could in principle excite molecules 
to a higher vibrational level, and a subsequent radiative 
cascade or vibration-vibration (V-V) transfer by near- 
resonant collisions could produce emission in the 3.25 fjxa 
bands. These processes are collectively known as non- 
resonant fluorescent pumping. An example is the fluores- 
cent e mission seen in the 10 /xm carbon dioxide bands on 
Mars ( Deming fc Mumma| [T983 ) . In that case, radiative 
absorption at short wavelengths is followed by rapid V-V 
collisional transfer to the upper state of the fluorescent 
emission. An analogous process might be energetically 
feasible for HD 189733b, since the stellar pumping could 
occur at short wavelengths where a higher stellar pho- 
ton flux is available, and the pumping could in principle 
involve multiple absorbing bands. 

A detailed line-by-line non-resonant fluorescence model 
for a realistic planetary atmosphere is beyond the scope 
of this paper, but we can evaluate the degree to which 
non-resonant fluorescence can augment the flux from res- 
onance fluorescence for the most relevant molecular can- 
didate (methane) . Non-resonant fluorescence of methane 
in planetary atmospheres ha s been discussed by|Drossart| 
et al.| (1999 [) (Jupiter \ and |Martm- Torres et al.| (|199^ ) 



(Earth); in particular, |Martin- Torres et al? "j "19980 evalu 
ate the 19 strongest vibrational bands of methane in our 



own atmosphere, demonstrating that the contribution 
to bands at 3.3 //m by radiative cascades from higher- 
energy states is insignificant. We performed a similar 
calculation by using the CH4@Titan database (2.6 mil- 
lion transitions) to identify all the methane bands with 
upper-state energies greater than 3 /mi that could then 
cascade downward; in total we found 188 bands. If we 
assume the absorbing bands are optically thin (appropri- 
ate for NLTE radiative pumping of higher-energy vibra- 
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tional bands), we can use the stellar flux at each wave- 
length from a Kurucz stellar model combined with the 
Einstein B coefficient of absorption for each transition 
to calculate the total amount of stellar flux that could 
be absorbed per second by gas at 1200K and compare 
it to the total flux absorbed by all the tabulated bands 
between 3 and 3.5 /im (a total of 52 bands). This cal- 
culation does not include the branching ratios for each 
upper state, which determine how the absorbed energy 
is split between specific lower states, and therefore pro- 
vides an upper limit to the fraction of absorbed energy 
that would actually cascade into the 3.3 /im region. We 
find that non-resonance pumping by stellar photons and 
subsequent radiative cascades could augment the effect 
of resonance fluorescence by a maximum factor of f .8. 
This contribution is clearly insufficient to explain the en- 
ergy difference between the impinging stellar radiation 
and the apparent emission detected by S10. We have 
similarly examined all our different molecular candidates 
scaled to the contrast claimed by S10 at 3.25 /im, and 
every one of them can be ruled out using the same com- 
putational approach. 

We therefore conclude that, even if the bright emission 
claimed by Sf is real, but variable so that it escapes our 
detection, it is too bright to be produced by fluorescent 
pumping from the stellar continuum. 

4.2. Time- Varying Emission from the Star or Planet 

Though the presence of gas-phase line emission can be 
definitively rejected in our data, the features detected 
by S10 could be attributed to a variable excitation or 
extinction mechanism that resulted in emission in 2007 
and a lack of emission in 2009. Regardless of the specific 
mechanism needed to generate the emission claimed by 
S10, we can also consider more generally whether or not 
changes in the properties of the star or the planet might 
account for our non-detection of the putative signal. 

One possibility that might produce highly variable 
emission is stellar flares, which could produce molecular 
excitation at high altitudes by charged particle interac- 
tions or enhanced UV pump ing. The star is active and 
variable at the percent level ( |Pont et al.|2007 1 , and solar 
flares have been known to produce variability in the auro- 
raeof giant planets in our own Solar System ( Pryor et al. 
2005 1 . We cannot rule out this possibility, but we be- 
lieve that any flaring event strong enough to produce the 
emission feature seen by SfO would produce enough high- 
energy radiation to dissociate the emitting molecules; ad- 
ditionally, there would be many other regions of the spec- 
trum that would show significant variability. A second 
scenario would be a decrease in the number of emitting 
molecules in the planet's atmosphere, possibly due to an 
increase in cloud opacity; however, no variability was de 



tected at 8 fim by |Agol et al] ( 2010 ) over seven secondary 
eclipses between 2006 and 2008 down to a level of 2.7%, 
suggesting that cloud variability over these timescales is 
minimal. 

A third source of variability which could be thought 
to mimic the effects of a planetary eclipse would be a 
change in the stellar radiation through the rotational 
modulation of stellar sp ots, which can a dd time-varying 
color-dependent effects ( Pont et al.|2008 ) , or a darkening 
of the planet's emission due to features on the planetary 
limb. However, we expect any modulations of the stellar 



or planetary brightness within a single night to be far 
below the level needed to produce the features seen by 
S10 (esp ecially given that the star's rotation period is 
11 days (Winn et al. 20071), and the variability would 
almost certainly not be restricted to a single spectral 
channel. 

4.3. Possible Broad Emission Features 

The absence of narrow emission lines in our high- 
resolution residual spectrum cannot rule out the presence 
of very broad emission features whose flux is spread over 
a large number of resolving elements in our spectra. The 
results presented by S10 are grouped into 0.1 /im spectral 
channels, while our entire spectrum covers only a frac- 
tion of one bin (~ 0.05 /im) at high resolution; therefore 
broad underlying features would be removed by contin- 
uum fitting in our reduction process. However, the pri- 
mary sources of broad spectral features - broadened gas- 
phase lines and PAH emission - can be rejected. Broaden- 
ing of individual molecular lines would need to approach 
130 km/s for this explanation to be viable, requiring ei- 
ther extreme pressure broadening or opacity b roadening. 
Based on broadening coefficients ( |Pine||1992[ ), about 20 
atmospheres of pressure would be needed; these high den- 
sities and pressures are incompatible with NLTE emis- 
sion, so pressure broadening can be excluded. Opac- 
ity broadening could result from extreme saturation of 
methane or water lines near 3.3 /im, similar to those 
commonly seen in transmittance spectra of Earth's at- 
mosphere (e.g., Fig. [5J near 3.280 //m), but would imply 
similarly bright emission at 2 /im and at 6 /im (water) or 
8 /im (metha ne) ; this emission has not been detec ted in 
earlier HST flSwain et al. | [2009| and Spitzer/IRS ( [Grill 
mair et al. 2008p observations, and the K-band results 
reported by S10 also show no such emission. Alternately, 
the emission feature could be intrinsically broad; in par- 
ticular, a broad emission band at 3.3 /im due to poly- 
cyclic aromatic hydrocarbons (PAHs) has bee n detected 
in a wide range o f astrophysical environments ( Allaman- 
dola et al. 119891). However, in this scenario PAH band 



emission would also be expected at longer wavelengths 
(6.2 /im, 7.7/im, 8.7/im, and 11.3/im), with similar or 
even brighter intensity than the 3.3 /im feature. These 
features are not detected in the Spitzer/IRS results. 

4.4. Telluric Contamination 

The final explanation we consider is that the excess flux 
reported by S10 may have been caused by systematic con- 
tamination due to imperfect removal of variable telluric 
absorption or emission. At L-band wavelengths, one of 
the largest sources of variability in flux over timescales 
longer than 15 minutes is the change in atmospheric con- 
ditions. Molecular species in the Earth's atmosphere 
affect observations in the L band in two ways: by re- 
ducing the stellar flux transmitted through the atmo- 
sphere (transmittance) and by emitting thermal radia- 
tion with a Boltzmann temperature characteristic of the 
lower atmospheric layers (radiance) . Radiance is an addi- 
tive component and must be subtracted from the stellar 
flux (usually accomplished by using an ABBA nodding 
sequence), while transmittance is corrected by dividing 
by a normalization factor based on a determination of 
the amount of flux lost (in our data reduction process, 
this is accomplished using models of telluric extinction). 
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Fig. 5. — A plot of the effects of changes in telluric water absorption during the observing night analyzed by |Swain et al.| | |2010[| (UT 
August f2, 2007). Telluric water vapor affects both the transmittance (how much star light reaches the detector] and the sky radiance 
(how much background thermal emission arrives from the atmosphere). A) Measurements of precipitable water vapor (PWV) using data 
from the tipping radiometer at the CSO on Mauna Kea (black), with the mean in-eclipse and out-of-eclipse values overplotted (red). B) 
We calculated the chang es in transmittance (A Trans.) and radiance (A Rad) for the measured change in water vapor, and compared 
them to the res ults fro m Swain et al. ( 2010J . While we cannot identify with certainty the cause of the differences between our data and 
the results from |Swain et"al!] | |^0iu| , the similarity between the effects of changes in sky radiance and the Swain et al. results is consistent 
with inaccurate removal ot variable telluric features. 



Several lines of evidence buttress the the hypothesis 
that the S10 result is due to inadequate removal of tel- 
luric contamination. First, the best-fitting molecular 
model in all L-band spectral channels presented by S10 
is an emission model of water at 300K (see Figure [3b 3 but 
we see no evidence of velocity-shifted emission lines m our 
data. While the characteristic temperature of the atmo- 
sphere of HD 189733b would be expected to be closer to 
1200K, the characteristic temperature for water in our 
own atmosphere is 296K. This suggests that a change 
in terrestrial water vapor during the planetary eclipse 
could be the cause of the observed variation in eclipse 
depth with wavelength. 

S10 analyzed data for a single eclipse observed on UT 
August 12, 2007 (M. Swain, private communication). 
Measurements of submillimeter opacity can be used as 



a prox y for precipitable water vapor (PWV; [Masson 
(1994)), and m easurements fro m the tipping radiometer 
on Mauna Kea ( Radford|2002 ) indicate that the average 
PWV was significantly higher (~ 9%) in the hours be- 



fore and after the eclipse, compared with the period of 
the eclipse (Figure [5j\.) . This variation would result in 
a higher telluric transmittance during eclipse (contribut- 
ing to a decrease in the uncorrected eclipse depth), but 
would also result in a lower telluric radiance (contribut- 
ing to an increase in the uncorrected eclipse depth). Us- 
ing the LBLRTM atmospheric modeling code, we have 
generated telluric transmittance and radiance spectra for 
two water vapor abundances - the average PWV during 
eclipse (1.18 mm) and the average PWV before and after 
eclipse (1.28 mm). To convert the radiance variation to 
a total flux we must multiply by the slit area; we use a 
square with dimensions of the slit width used by S10 (1.6 
arcseconds). In Figure[5j3, we show the change in the tel- 
luric radiance emission m units of the stellar flux (ARad) 
and the change in stellar flux due to extinction by the tel- 
luric transmittance function (ATran). The fluxes shown 
for each spectral channel are the contributions from each 
effect on the apparent eclipse depth in that channel, if 
not compensated for in the original analysis by S10. 
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The similarity between the signal observed by S10 and 
the contribution from the sky radiance due to changes in 
the telluric water vapor emission suggests that the effects 
of atmospheric transmittance may have been removed ac- 
curately (possibly by normalizing each spectrum by the 
mean baseline flux) but that the effects of telluric ra- 
diance may have been incompletely corrected, possibly 
because S10 processed their A and B beams separately. 
Incorrect compensation for the effects of sky radiance 
would lead to an apparent drop in flux during the pe- 
riod of the eclipse and result in a spurious measurement 
for the loss of flux from the planet. This effect may 
have been exacerbated by the wide slit width used, which 
would have made accurate removal of the sky emission 
more difficult. The apparent ability of the authors to 
match the previous HST observations at K-band wave- 
lengths further supports radiant sky emission as opposed 
to a problem with transmittance correction, since the 
intensity of telluric radiance at K-band wavelengths is 
much lower than in the L band. 

Eventually it should be possible to directly assess the 
contribution of telluric contamination in the S10 results, 
but that assessment will require re-analysis of the S10 
data using their original methodology, which uses non- 
standard analysis techniques that are described in insuffi- 
cient detail in their paper to be reproduced, and compar- 
ison of these results with an analysis using independent 
methodology. Although our observations and analysis 
unequivocally contradict the S10 results, our analysis al- 
gorithms - being optimized for high spectral resolution 
data - do not allow us to test the specific methodology 
used by S10 with equal rigor. While we cannot confirm 
sky emission as the source of the feature reported by S10, 
in the absence of a viable astrophysical explanation we 
consider it as the most likely possibility. 

5. CONCLUSIONS 

We report a non-detection of emission at 3.3 /im from 
the extrasolar planet HD 189733b using high-resolution 
spectroscopic observations taken with the NIRSPEC in- 
strument on the Keck-II telescope, a result in contrast 
with the detection o f strong emission a t these wave- 
lengths announced by |Swain et al.| ( |2010[ ). We modeled 



the expected signal using ro-vibrational emission mod- 
els of molecular species with transitions in the relevant 
wavelength region using a wide range of rotational ex- 
citation temperatures and vibra tional level intensit ies, 
scaled to the results presented by Swain et al. (2010), in 
order to predict the line flux required at high resolution. 
No flux was present at any of the expected transition fre- 
quencies between 3.27 and 3.31 /im, with upper limits 40 
times smaller than the expected line fluxes. The condi- 
tions that would lead to broad emission features beyond 
our detection limits are extremely difficult to reconcile 
with realistic models and previous observations. Addi- 
tionally, our analysis indicates that the emission, if real, 
is too bright to be produced by fluoresence. Our wave- 
length region only covers a small section of the spectrum 
published by S10, and we cannot rule out an emission 
mechanism that produces flux outside our band pass. 
Additionally, we cannot rule out an exotic highly time- 
variable stellar emission process such as charged-particle 
excitation due to flares. However, we regard these ex- 
planations as improbable, and conclude that inadequate 
telluric correction is the most likely explanation for the 
Swain et al.l (I2010I) results. 
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